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Abstract 

We present a criterion of reducibility of a zero curvature representation 
to a solvable subalgebra, hence to a chain of conservation laws. Namely, 
we show that reducibility is equivalent to the existence of a section of the 
^0 ' generalized Riccati covering. Results are applied to conversion between 

■ Guthrie's and Olver's form of recursion operators. 

• i— I ' 

1 Introduction 

To establish integrability of a nonlinear partial differential equation in the sense 
of soliton theory [1, 37], at least in two dimensions, one usually looks for a zero 
\ curvature representation (ZCR) [43], possibly in the form of a Lax pair [17]. 

C^) • If depending on a non-removable (spectral) parameter, a ZCR may serve as a 

\& , starting point of methods to derive infinitely many independent conservation 

laws and large classes of exact solutions. 

However, certain ZCR's do not imply integrability because of specific degen- 
eracy, which does not even rule out possible dependence on one or more non- 
removable parameters. E.g., Calogero and Nucci [3] gave a formula to assign a 
Lax pair to any nonlinear system possessing a single conservation law, arguing 
that such systems are too abundant to be all integrable. Recently Sakovich [33] 
observed that the Calogero-Nucci examples can be singled out by properties of 
their associated cyclic bases. In particular, the 'bad' ZCR's fail to generate an 
; "*] ■ integrable hierarchy. 

In this paper we postulate that a ZCR is degenerate if it takes values in a 
solvable Lie algebra or is gauge equivalent to such. Even though some researchers 
are inclined to admit the relevance of such ZCR's to integrability, our results 
below seem to support the opposite opinion. 

ZCR's taking values in an abelian algebra are well known to be equivalent to 
a set of local conservation laws (see [1, Sect. 3.2.c]). Using the Lie theorem on 
finite-dimensional representations of solvable algebras, we show in Sect. 4 rather 
easily that every ZCR that takes values in a solvable algebra is equivalent to a 
'chain' of nonlocal conservation laws. This simple result renders, e.g., attempts 
to generate infinitely many independent conservation laws out of a degenerate 
ZCR rather unrealistic. 

In Sect. 5 we address the problem of detecting reducibility of a ZCR to a 
subalgebra, in particular, to a solvable one. Purely algebraic criteria are insuf- 
ficient since the Lie algebra a ZCR takes values in may be altered by gauge 
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transformation. On the other hand, when trying to find the reducing gauge ma- 
trix directly, we face a rather large underdetermined differential system. Our 
idea is to employ an appropriate matrix decomposition, namely, the Gram or 
Gauss decomposition. Earlier these decompositions were applied by Dodd and 
Paul [6, 7] in the context of Backlund transformations. A remarkable connection 
between decompositions and integrable systems emerged in numerical analysis 



The last section is devoted to recursion operators, direct and inverse, for sym- 
metries of integrable systems [2, 14, 28, 39]. In Olver's [27, 28] formalism, a 
recursion operator is a linear integro-diffcrential operator 'J, which maps sym- 
metries to symmetries. The standard way of inverting \1/ consists in finding dif- 
ferential operators K,L such that \& = LoK^ 1 ; then "f -1 = K oL^ 1 . However, 
one encounters the problem of writing the inverse L^ 1 as an integro-differential 
operator. In the scalar case, L may be put in the form L = q n Dq n -\D ■ ■ ■ qiDqa, 
where the coefficients qt are expressible as quotients of wronskians of indepen- 
dent solutions Vi of L(v) = (see [31] for a simple derivation of this classical 
result, equivalent to decomposability into first-order factors; see [44] for the 
matrix case). In our context, are nonlocal functions and finding them is con- 
sidered to be the most difficult part of the whole procedure. Once q^ are found, 
one can invert L simply as = q^ ) 1 D~ 1 q^ 1 ■ ■ ■ D~ 1 q~^ l D~ 1 q~ 1 . This is es- 
sentially the general scheme behind the works [12, 18, 19, 20, 29, 32]. 

Guthrie's [11] recursion operators resemble Backlund autotransformations for 
the linearized system and indeed can be interpreted this way (see [22]); their 
inversion is quite straightforward and does not require the introduction of new 
nonlocalitics. Moreover, Guthrie's operators do not suffer from the known ab- 
normities, related to the fact that D~ 1 oD = id fails to hold ([10, 34]). Let us also 
remind the reader that computing the 'inverse' Guthrie operator starting from a 
known ZCR may turn out to be easier than computing the 'direct' one (see [24]). 
The conversion from Olver's to Guthrie's form was explained by Guthrie [11] 
himself, the result being further strengthened by Sergyeyev [35]. Concerning 
the backward conversion, the x-part of a Guthrie operator can be written as an 
integro-diffcrential operator if the ZCR underlying it is lower triangular. A non- 
parametric ZCR can be made lower triangular at the cost of the introduction 
of appropriate nonlocalities. To introduce only few (respected) nonlocalities, we 
take into account a particular observation (already exploited in [24]) about the 
structure of Guthrie's recursion operators of integrable systems. 

2 Preliminaries 

Let E be a system of nonlinear partial differential equations (PDE) 



on a number of functions u k in two independent variables x,y. Here each F l is 
a smooth function depending on a finite number of variables x, y, u k ,u k , u k , . . ., 
Uj,. . . , where I stands for a symmetric multiindex over the two-element set of 
indices {x,y}. Besides the local variables x,y,u k ,u k , we shall also need non- 
local variables or pseudopotentials [40], which may be introduced as additional 
variables satisfying a system of equations 



[4, 5, 41]. 



F l = 
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where /*, g % are functions depending on a finite number of local variables as well 
as the pseudopotentials z 3 ; we require that the system (2) be compatible as a 
consequence of (f ). 

Within their geometric theory of systems of PDE's, Krasil'shchik and Vino- 
gradov [15] introduced the notion of a covering, which separates the invariant 
content of nonlocality from its coordinate presentation. Pseudopotentials then 
correspond to a particular but arbitrary choice of coordinates along the fibres of 
the covering in question. We recall the basic facts below; details we had to leave 
aside may be found in [15] and also in [16, Ch. 6]. Readers interested mainly in 
practical computations may skip the rest of this section. 

Let J°° be an infinite jet space equipped with local jet coordinates x, y, u , Uj, 
the functions F l then may be interpreted as functions defined on J°°. Since all 
our considerations are local, we simply let J°° be the space of jets of sections of 
the trivial fibred manifold Y x M — ► M, where M = M. 2 , with x, y being coordi- 
nates on R 2 and u k being coordinates on Y. On J°°, we have two distinguished 
commuting vector fields 



which are called total derivatives. 

The equation manifold £ associated with system (1) is defined to be the sub- 
manifold in J°° determined by the infinite system of equations F l = and 
D]F l — for I running through all symmetric multiindiccs in x,y. The total 
derivatives D x ,D y are tangent to £, therefore they admit a restriction to £. 
In what follows, equations will be identified with equation manifolds equipped 
with the restricted total derivatives; this approach is indeed very practical and 
suitable for all needs to be encountered below. 

Mappings between equation manifolds that commute with projections to the 
base manifold M and preserve the total derivatives will be called morphisms of 
equations; they map solutions to solutions (we shall not use the general mor- 
phisms of diffieties which need not commute with the projections and only pre- 
serve the distributions generated by the total derivatives). Bijective morphisms 
are called isomorphisms; their inverses are isomorphisms, too. 

A covering over an equation £ consists of another equation £' and a surjective 
morphism £' — ► £. 

The system formed by Equation (1) and the 2k additional equations (2) gen- 
erates a covering, where £' is the trivial vector bundle £ x K fe and z 1 , . . . , z 
provide coordinates along K fe . In particular, the projection preserves the coor- 
dinates x, y. If /*, g l are functions defined on E' such that the vector fields 



commute (which is a geometric way of saying that Equations (2) are compatible), 
then £' equipped with the vector fields (3) is a /c-dimensional covering over £. 
Recall from [15] that every finite-dimensional covering is locally of this form. 

Two coverings £' and £" are said to be isomorphic over £ if there exists 
an isomorphism of the equations £' = £" that commutes with the projections 
to £. Isomorphic coverings result from invertible transformations of nonlocal 
variables. A fc-dimensional covering is said to be trivial if it is isomorphic to one 
with p — g % — 0; such a covering is essentially a family of identical copies of £. 





(3) 
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The simplest yet useful covering (2) may be associated with a single nontrivial 
conservation law a = f dx + g dy, i.e., a pair of functions /, g defined on £ and 
satisfying D v f = D x g on £: 

Definition 1 A one- dimensional abelian covering associated with a conserva- 
tion law a = f dx + g dy is defined to be the trivial vector bundle £ x M — ► £, 
equipped with total derivatives 

D' x = D x + fl z , D' y = Dy+4 z , 
where z denotes the coordinate along R. 

As /, g do not depend on z, the vector fields D' x ,D' y on £' commute if and 
only if D y f = D x g. The variable z is called the potential of the conservation 
law a. We have D' x z = f, D' y z = g or briefly z x = f, z y = g. 

Recall that a conservation law is said to be trivial if there exists a (local) 
function h on £ such that / = D x h, g = D y h. A covering associated to a trivial 
conservation law is isomorphic to a trivial covering through the invertible change 
of variables z — z' + h. 

A covering £ — ► £ is said to be trivializing for a conservation law a = f dx + 
g dy, if the pullback a of a along the projection £ — > £ is a trivial conservation 
law on £ . Obviously, the one-dimensional abelian covering associated with the 
conservation law a trivializes a. 

A general n-dimcnsional abelian covering is obtained by repeating the con- 
struction of the one-dimensional abelian covering (cf. [40, Sect. IV]): 

Definition 2 An n-dimcnsional covering £ over £ is said to be abelian, if 

(1) either n = 1 and £ is a one-dimensional abelian covering over £ in the 
sense of Definition 1; 

(2) or £ is a one-dimensional abelian covering over an (n — l)-dimcnsional 
abelian covering £' over £. 

Let us remark that Khorkova [13] introduced the universal abelian covering, 
which need not be finite-dimensional. 



3 Zero-curvature representations 

Simplest pscudopotentials that are not potentials are associated with non- 
degenerate zero-curvature representations. Let g be a matrix Lie algebra (recall 
that according to the Ado theorem every finite-dimensional Lie algebra has a 
matrix representation). By a Q-valued zero- curvature representation (ZCR) for 
£ we mean a q- valued one-form a — A dx + B dt defined on £ such that 

D y A - D X B + [A, B] = (4) 

holds on £, which means that (4) holds as a consequence of system (1) (we do 
not insist that (4) necessarily reproduces system (1), which is normally required 
in integrability theory). 

Let Q be the connected and simply connected matrix Lie group associated 
with q. Then for an arbitrary (^-valued function S, the form a s = A s dx + B s dt, 
where 

A s = D X SS- 1 + SAS- 1 , B s = DySS- 1 + SBS- 1 (5) 
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is another ZCR, which is said to be gauge equivalent to the former. 

A ZCR is said to be trivial if it is gauge equivalent to zero, i.e., if A = DxSS^ 1 , 
B = DySS^ 1 . A covering £ — ► £ is said to trivialize a ZCR a = Adx + B dy 
if the pullback a of a along the projection £ — ► £ is a trivial ZCR. 

A trivializing covering for the ZCR a can be obtained in the following way. 

Proposition 3 For every g-valued ZCR a on £ there exists a covering ir a : 
£ a — ► £ that trivializes a. 

Proof Let a = A dx + B dy be a ZCR, where A and B are n x n matrices 
belonging to the algebra g. Put £ a = £ x G, where G is the matrix Lie group 
associated with g. Given an element C G fl, let us denote by £c the right- 
invariant vector field on G corresponding to C. Given a g- valued function C on 
£, let us denote by He the unique vector field on £ a with the ^-component zero 
and the G-component equal to £c, at each point of £ a . Considering the vector 
fields 

D x = D x + Ha, D y = D y + Hb 

on f a , where D x ,D y are the total derivatives on £, let us show that D x ,D y 
are the total derivatives for a trivializing covering ir a : £ a — ► £ of a. 

Let A = (dij), B = (bij). Let us first consider G = GL„ with its natural 
parametrization GL„ = {(zij) | dctz^ ^ 0}. We have 

Ha = ^2 aij Zji — , H B = ^ b l0 Zji~^-- 

Then D x ,D y commute since 

[D x , D y ] = [D x ,D y ] + [D X ,E B ] - [Ha, D y ] + [Ha, H b ] 

= ^D x B — D y A — [A,B] 

= 0. 

The same holds for arbitrary G C GL n , since the vector fields Ha , c b are 
tangent to G whenever A,B belong to g. 

Now denote by 6 the projection £ a = £ x G — ► G viewed as a matrix-valued 
function on £ a . Then D X Q = and therefore 

(D x O)^ = (H A e) M „ = ^aijZji- — z^ v = ^2,a^z jv = (AQ)^. 

, &Zn 

Thus, D X Q ■ 6 _1 = A and similarly D y Q ■ Q^ 1 = B, whence the pullback of a 
on £ a is trivial. 

The system (2) corresponding to £ a can be compactly written in terms of a 
single matrix 9 as 

9^ = AQ, O y = BS. (6) 

Under the gauge transformation (5), the matrix 6 becomes SQ. The coverings 
£ a and £ a s are isomorphic via i — ► SQ. 

The trivializing covering ir a just constructed has the following factorization 
property: 
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Proposition 4 Let p : £' — ► £ over M be a trivializing covering for a ZCR a 
on £ . Then there exists a morphism p" : £' — ► £ a such that n a o p$ = p. 

Proof Let a = Adx + Bdy. Since p is over M, we have p*a — p*Adx + 
p*B dy. By assumption this is a trivial ZCR, whence p*A = D^SS" 1 and p* B = 
D'ySS^ 1 for a suitable G-valued function S on £'. Recall that fibres of the 
covering £ a arc diffcomorphic to the Lie group G. Therefore we can define 
a mapping p" : £' — > £ a by the formula 6 o p" = 5, where, as above, 6 
denotes the projection £ a = £ x G — ► G. The mapping p" is a morphism, since 
(eop«) x = 5^ = ,45 = A6op«. 



4 Lower triangular ZCR's 

Let in denote the algebra of matrices 

\ 



/an 

0-21 0-22 

031 «32 033 

Wnl a n2 a«3 
(fc) 







(7) 



Denote by t„ , fc > 1 , the derived algebra formed by matrices satisfying ay = 
whenever i — j < fc. 
ZCR's with values in t„ are, in a sense, equivalent to an abelian covering. 

Proposition 5 Every i n -valued ZCR can be trivialized by means of an abelian 
covering of dimension < \n(n + 1). 

Proof Let a = A dx + B dy be a ZCR such that matrices A and B are lower 
triangular. We shall construct an abelian covering £( n_1 ) in n steps. 

It follows from Equation (4) that 71 = an dx+bu dy, 72 = a 2 2 dx+b 2 2 dy, . . . , 
In = a nn dx + b nn dy are conservation laws. Let us denote by £^ the associated 
abelian covering with potentials h\,...,h n satisfying 



1, . . . ,n. 









bu 


for i = 


1, 


Then 

























\ 







e -h 2 










H = 








e -hs 









V 








. e -h, 


/ 



is a matrix defined on £ (°\ with the property that all diagonal entries of the 
gauge equivalent matrix A' = A H vanish: 



A' 



( 

2 31 
\ a 'nl 



l 32 



l n2 



n,n — l 







(8) 
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and similarly for B'. Hence, A' , B' take values in tnK 

By the same Equation (4), 7 2 = a' 2l dx + b' 2l dy, -f' 3 — a' 32 dx + b' 32 dy, . . . , 
In = a ' n _ ln dx + b' n _ ln dy are conservation laws on £^°\ Let us introduce a 
covering £' over £^> with potentials h' 2 ,...,h' n satisfying 

h'i,x = a i,i-i> Ky = b 'hi-i fori = 2,...,n. 
Denoting 



/ 1 



H' = 



-h'o 



-ho 1 



V 











-tin 1/ 



we see that the gauge equivalent matrices A" = A ,H and £?" = B ,H take 
values in t n ' now. Compared with (8), A" and £r have one more subdiagonal 
of zeroes. The next step is similar: 73 = a' 3 \ dx + b 31 dy, 7" = a![ 2 dx + b'{ 2 dy, 



dx 



b n -2,n 



dy are conservation laws on £' . Let us introduce a 



covering £" over £' with potentials h 2 , . . . , /i" satisfying 



h'J 



-2' 



for 



, n. 



Denoting 



H" = 



( 1 











1 

1 

-til 



V 



we observe that A" 



\nH" 



0\ 





B" 



-< l) 



B 



tlH" 



take values in tn\ and so on. 



Continuing the process until A^ , B^ become zero, we end up with a sequence 
of \n{n + 1) conservation laws 



7i 



72 

72 



73 
73 

7 3 ' 



(ra-2) 
7n-l 



7n 
7; 
7* 



(9) 



7° 

In 

7, 



2) 

(n-1) 



where (a) 71, . . . , 7„ are conservation laws on £; (b) 7^ t ^ 1; . . . , 7!" ^ are con- 
servation laws defined on the abelian covering £(™ _t_1 ) associated with the 
conservation laws of all the previous levels. 

Finally, a HH " H ' = = 0, where each ifM is defined on £M. Summing 
up, the covering trivializes a. 



The sequence (9) will be called an n-fold chain of conservation laws. 
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Proposition 6 Let a be a t n -valued ZCR, then the associated covering ir a is 
isomorphic to an abelian covering of dimension < ^n(n + 1). 

Proof According to Proposition 5, there is an abelian coverings : — ► £ 
that is trivializing for a; namely, we have a K = 0, where K = HH' ■ ■ ■ i^™ -1 ) 
(see proof of Proposition 5). Hence, a = K and, according to Proposition 4, 
there is a morphism : ff™ -1 ) — > £ a , given by 8 = K~ x . Here repre- 
sents the totality of coordinates along the fibres of the covering £ a , while K 
is parametrised by coordinates along the fibres of the covering £( n ~ 1 \ It 
follows that p" is bijective on the fibres, hence isomorphism. 

5 Reducibility 

A g- valued ZCR is said to be reducible if it is gauge equivalent to a ZCR taking 
values in a proper subalgebra f) C q; otherwise it is said to be irreducible. 

Let f) C g be a subalgebra. We present a simple criterion for reducibility of 
a g- valued ZCR to t). Let H c Q be the Lie subgroup corresponding to the 
subalgebra f). We call H a right factor if there exists a submanifold JC C G 
(possibly with singularities) such that the multiplication map 

fi-.ICxH^G, (K, H) i — ► KH (10) 

is a surjective local diffeomorphism. The manifold JC will be called a cofactor. 
By surjectivity, every element S £ Q can be decomposed as a product S = KH, 
where K e JC and H e H, possibly non-uniquely. The map fi being a local 
diffeomorphism, JC has the minimal possible dimension dim/C = dimQ — dim7i. 
If H is closed, then the assignment K i — ► KH defines a local diffeomorphism of 
JC onto the homogeneous space Q/H. 

Proposition 7 Under the above notation, a Q-valued ZCR a on £ is reducible 
to the subalgebra f) if and only if there exists a local JC-valued matrix function 
K on £ such that a K lies in f). 

Proof The gauge equivalence with respect to H £ H preserves the subalgebra 
f). Therefore, the gauge-equivalent ZCR a s — (a K ) H lies in f) if and only if a K 
lies in f). 

Otherwise said, if a ZCR is reducible to (), then the corresponding gauge 
matrix can be found in JC. Understandably, different choices of the cofactor JC 
may lead to different reducibility criteria. 

In this paper we are primarily interested in detecting reducibility to a solv- 
able subalgebra. By the well-known Lie theorem ([9, Sect. 9.2]), every finite- 
dimensional representation of a solvable Lie algebra is equivalent to a represen- 
tation by lower triangular matrices. Hence, every ZCR reducible to a solvable 
subalgebra is reducible to t„ (and can be trivialized using an abelian covering 
according to Proposition 5). 

Let us therefore apply Proposition 7 to \) = t n . There are two standard ways 
to make t„ a right factor in gl n . 

The QR or Gram decomposition is an alternative formulation of the fa- 
mous Gram-Schmidt ortogonalization algorithm. Namely, every n x n complex 
matrix A can be decomposed as a product A = QR, where Q € SU„ and R G t„ 
[25, Ch. 6, Sect. 1.9]. Proposition 7 then yields 
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Proposition 8 A real (complex) ZCR a on £ is reducible to lower triangular 
if and only if there exists an SO n -valued (SU n -valued) function K on £ such 
that a K is lower triangular. 

However, the factors Q and R are unique up to a unimodular diagonal mul- 
tiplier: QR = QQ ■ Q^R, where 6 = diag(0i, . . . , 6 n ) e S(Ui x ••• x Ui), 
i.e., \9 L \ = 1 and IIILi &t = 1- Thus, the mapping (10) is not a local diffeo- 
morphism unless it is restricted to a suitable immersion of the orbit space 
SU„/S(Ui x • • • x Ui) into SU„. In the real case we have 6 L = ±1 and we 
get a 2"~ 1 -to-one local diffeomorphism (10) with JC = SO„. 



The LU or Gauss decomposition can be derived from the Gaussian elim- 
ination algorithm. The following result is well known ([25, Ch. 6, Sect. 1.8]): 

Proposition 9 For every non-singular matrix A there exist matrices P, U, L 
such that A = PU L, L is lower triangular, U is upper triangular with diagonal 
entries equal to 1, and P is a permutation matrix. The matrix P can be omitted 
if and only if all principal minors of the matrix A are nonzero. 

(Recall that Gaussian elimination may require row swapping, which is where 
the permutation matrix P comes from.) Let K, denote the set of all products 
PU where P is a permutation matrix and U is an upper triangular matrix with 
diagonal entries equal to 1. Then if is a union of n! intersecting submanifolds, 
labelled by permutation matrices P. Compared to the QR-decomposition, each 
of the n! submanifolds is easier to parametrize than SO or SU. 

Proposition 10 A ZCR a on £ is reducible to lower triangular if and only if 
there exists a permutation matrix P and a matrix-valued function 



H 



/I hi2 h 13 
1 h 23 
1 



(11) 



Vo • -oi/ 

on £ such that a PH is lower triangular. 

Before looking more closely at low values of n, we make a general remark to 
the effect that every g[„-valued ZCR is reducible to st n : 

Remark 11 A g[„-valued ZCR is decomposable into an s[„-valucd ZCR (trace- 
less summand) and a conservation law (the trace). 



5.1 The case of n = 2 

When n = 2, the reducibility condition corresponding to the QR decomposition 
is: 



Proposition 12 A Qi 2 -valued ZCR 
'an ai2 



a = Adx + B dy 



«21 022 



dx 



hi h 2 

&21 b22 



dy 
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is reducible to lower triangular if and only if there exists a function 4> on E that 
is a solution of the system 

D x <p = —a\2 cos 2 (f> + (an — 0,22) sin cos + (Z21 sin 2 <p, 

(12) 

D y (f> = —b\2 cos 2 + (6n — b 22 ) sin cos <j> + &21 sin 2 0. 



Proof An arbitrary SO2 matrix is 

/ cosd> sin< 
\—sm<p cost 

By Proposition 8, the ZCR a is reducible to lower triangular if and only if a K 
is lower triangular, which is exactly the meaning of conditions (12). 

The reducibility conditions corresponding to the LU decomposition are: 
Proposition 13 A Qi 2 -valued ZCR 

a = Adx + B dy = ( ° n CLl2 \dx+\ 1 11 ^ 12 )dy 
\a-21 0,22 J \02l o 22 / ' 

on E is reducible to lower triangular if and only if 
1. either there exists a local function p on £ such that 



D x p = -ai2 + (an - a 22 )p + a 2 ip 2 , 
D y p = -612 + (fen - b 22 )p + &21P 2 ; 
2. or A, B are upper triangular: 

021 = 621 = 0. 

Proof An arbitrary /C-valued function is K = PU, where 



(13) 



U 



1 p 
1 



and P is one of the two 2x2 permutation matrices 

Pl2 = (o 1)' P21 = (i 

Subcases 1 and 2 correspond to the choices P = Pi 2 and P = P 2 i, respectively, 
and express the conditions that A PU , B PU be lower triangular. 

Recall that a quadratic or Riccati pseudopotential p associated to an g[ 2 -valued 
ZCR a is defined by the compatible system 

Px = -012 + (an - a 22 )p + a 21 p 2 , 

(14) 

P y = -&12 + (&11 - b 22 )p + b 2 ip 2 , 

which are essentially Equations (13). The system (14) being compatible, let us 
introduce the corresponding one-dimensional Riccati covering. Proposition 13 
then says that a non-upper-triangular ZCR is reducible to lower triangular if 
and only if the Riccati covering has a local section. 
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Remark 14 Obviously, Equations (13) and (12) are not independent, the ex- 
plicit mapping of their solutions being p = tan0 for 4> 7^ (2fc + l)7r/2. Recently 
Reyes [30] pointed out a geometric interpretation of the same correspondence 
in terms of pseudospherical equations. 

Example 15 The Burgers equation u t = u xx + uu x is well known to be inte- 
grable via the Cole-Hopf transformation, which relates its solutions with those 
of the heat equation [1, Sect. 3.1]. Of the several Lax pairs that have been found 
all turn out to be degenerate. Let us consider one example [26, 42], where the 
lower triangular representation could not be obtained by purely Lie algebraic 
methods: 



a 



1 

- \ u x + i ( u + A) 2 



dx 



+ 



4 ^xx 



\{u- \)u x + ^ (u - A)(u + A) 

In this case, Equations (13) have a local solution p = 4/(u + A), hence 

1 4/(« + A) 
1 

is a gauge matrix to make the ZCR a lower triangular. 

That the Burgers equation has no irreducible g[ 2 -valued ZCR follows from 
the recent classification of second-order evolution equations possessing an sfo- 
valued ZCR [23] and Remark 11. The non-existence of an irreducible ZCR of 
Wahlquist-Estabrook type for arbitrary n was proved Dodd and Fordy [8] who 
established solvability of the Wahlquist-Estabrook prolongation algebra of the 
Burgers (and also of the Kaup) equation. 

Example 16 The Calogero-Nucci example [3] of a ZCR that exists for every 
equation possessing a conservation law f t = g x : 







fx 



+ 



VT + A/5 



■W9- 



2 9 9x 

v l 7 



dx 



9_ 

f 



V9 



"1 



(15) 



dy 
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where 77, A, /1, v are arbitrary constants. This ZCR is reducible, which follows 
from Proposition 13 along with explicit formulas for its reduction. Indeed, we 
have Subcase 1 again and one easily finds a local solution 



= 1 (m+ vV + 4A),/--2t7 
P ~ 2 Xf + W f - n 2 

of Equations (13). Hence, the above ZCR is reducible to lower triangular. 
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Continuing the reduction further, one finally arrives at an abelian subalgebra. 
Namely, if p is as above and 

= A/ 2 + wf - V 2 



Vm 2 + 4A / 



dx 



(A/ 2 + W f - 77 2 )(2AJ + vV + 4A)»?) 

v = ^ 

2A/ + (/i+ vV + 4A)?7 

then the product of gauge matrices 

(v^// \/l 0\/l p 
1/V?)\q l)\0 1 

takes the ZCR to the diagonal form 

A(M- vV+4A)/-r/ J) ^ 

V V / M 2 + 4A)/-? ?y 

+ Vm 2 + 4A) 5 + z/ \ 

V v 7 ^ 2 +4A).g + i// 

which is manifestly equivalent to the conservation law f dx + g dt. 
5.2 The case of n > 3 

For n > 3, the QR approach is impractical due to relative clumsiness of the 
parametrisation of SO„ by generalized Euler angles. On the other hand, the LU 
criteria come out subdivided into as much as n\ subcases, one for each of the n\ 
permutation matrices P. 

For every n, the case of general position occurs when all principal minors 
of the gauge matrix K are nonzero. Then the permutation matrix P equals 
the identity matrix and we can derive explicit formulas that generalize (13) to 
arbitrary n. 

Proposition 17 A gl n -valued ZCR a = Adx + Bdy, where A = (ay) and 
B = (bij), is reducible to lower triangular by means of a gauge matrix with 
nonzero principal minors if and only if the system 



D x pki= - (- 1 Y a ki B Pt„i 1 Pi 1 i 2 ---Pt 



0<r<n-l 
io<ii<---<V=i 



^ ' ( l) Pkj a jioPioiiPiii 



' Pi r -\i r I 



0<r<n-l 
k<j 
~i<il<---<i r =l 



D y Pkl= ~ ^2 (- 1 ) r bki Pi iiPiii2 ■ ■ -Pir-iir 



(16) 



0<r<n-l 
io<il<-<i r =l 



^ ' ( 1) PkjbjioPiohPiih ' ' ' Pi, — i*r 

0<r<n-l 
k<j 

ia<il<---<i r =l 

on \{n~ l)n unknown functions pu, k < I, has a local solution. 
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Proof According to Proposition 9, every gauge matrix S with nonzero principal 
minors decomposes as S = LU, with L lower triangular and 



U = 



/l Pl2 P\3 
1 P23 

1 



\0 

The inverse of U is 



Pln\ 

P2n 
Pin 

1 J 





(I 


912 


913 


• • • qin\ 







1 


923 


■ ■■ q_in 










1 


■ ■ ■ 9371 




V° 








... l) 


where 










Qij = 




E 


(- 


i) r Pi iiPi 



l<r<n-l 

t=io<*l<'"<»r = 7 



(-1) 



i+j 



Pi,i+1 Pi,i+2 

1 Pi+l,»+2 Pi+l,»+3 
1 Pi+2,i+3 Pi+2,i+4 

1 







Pi,j-1 



Pi,j 
Pi+l,j 



Pj-2,j-l Pj-2,j 
1 Pj-lJ 



since + ^^KiKlPkiQu + Pki = whenever k < I. Let us consider the gauge 
equivalent matrix A f ' : = U x U~ l + UAU~ l . Terms that contain total derivatives 
.D^py can occur only in the first summand, which is 



Z23 











\0 



Zln\ 

Z2n 



where 

Zkl 



E 



(-l) r 1 D x p ioil -p ili2 ■ ■■p lr _ ll 



l<r<n-l 
k=»o<*i<"-<ir=J 



= (-D 



k+l + l 



D x Pk,k+l D x pk^+2 ■ D x pk,l-1 DxPk,l 

1 Pk+l.k+2 Pk+l,k+3 • Pk+l,l 

1 Pk+2,k+3 







Pl-l,l 
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for all k < I. Denoting A u =: A' = (a^-), we have 

a 'ki : = z ki + an + ^2 akjqji + *^2pkidijqji + ^Pfcjflti- 

The condition of A' being lower triangular, a' kl — for all k < I, constitutes a 
system of equations in total derivatives D x pij. The equivalent system resolved 
with respect to the derivatives is a' k[ +^2 k<h<l a' kh Phi — 0, since derivatives occur 
only in the summands containing Zij, which are z ki + J2k<h<i z khPhi = D x p H . 
The remaining summands then simplify to the expressions given in the statement 
of the proposition. 

5.3 The generalized Riccati covering 

A tedious computation shows that Equations (16) are compatible, meaning that 
there are no integrability conditions resulting from the equalities D xy pki = 
DyxPki- This implies the existence of a covering associated with a ZCR which 
naturally generalizes the Riccati covering. Similar result holds for more general 
types of decomposition, too. 

Let a subalgebra 6 C g be a direct complement to the subalgebra t) C g 
considered throughout this section. Let pr { : q — > t be the corresponding pro- 
jection. Then the condition a K e t) (Proposition 7) can be equivalently written 
as pi t a K = 0. 

Denoting by JC the Lie connected and simply connected matrix Lie group 
associated with the subalgebra 6, we have 

Proposition 18 Under the above notation, the differential equations 

Wt {K x K- 1 +KAK- 1 ) = 0, 

WtiKyR- 1 +KBK- 1 ) = 
on a matrix K e JC are compatible. 

Proof Since K S JC, matrices KxK" 1 , KyK^ 1 belong to 6 and are mapped iden- 
tically under the projection pr { . Hence Equations (17) are differential equations 
on K and, moreover, can be resolved with respect to K x ,K y . Let us consider 
their derivatives 

= By Wl {K x K- x + KAK- 1 ) = Wl {K xy K- x - KBAK^ 1 + KAyK^ 1 ), 
= D^iKyK- 1 + KBK- 1 ) = ^{Ky^- 1 - KABK' 1 + KB^ 1 ), 

where we have made substitutions pr^KxK^ 1 ~> —pr^KAK^ 1 , pr t K y K~ 1 ~» 
— px^KBK^ 1 according to (17). These equations can also be resolved with re- 
spect to K xy and K yx , respectively. Now one can perform the standard check 
that K xy coincides with K yx : 

Wt(K xy - K yx )K- x = pi t K{A v ~B X +AB- BA)^ 1 = 0. 

Definition 19 Given a ZCR a of an equation £ and the decomposition q = h+J 
as above, we define the associated generalized Riccati covering as £ x JC — ► £, 
assuming that the corresponding matrix of pseudopotentials K e JC satisfies 
Equations (17). 
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Summing up, we obtain: 

Corollary 20 A gl n -valued ZCR a is reducible to lower triangular by means of 
a gauge matrix from IC if and only if the generalized Riccati covering associated 
with the decomposition g[„ = t n + 6 has a local section. 

Choosing t to be the Lie algebra of strictly upper triangular matrices, we have: 

Corollary 21 A gl n -valued ZCR a is reducible to lower triangular by means of 
a gauge matrix with nonzero principal minors if and only if there exists a local 
solution to Equations (16). 



6 Guthrie's formulation of recursion operators 

In 1994, G.A. Guthrie [11] suggested a general definition of a recursion op- 
erator, free of some weaknesses of the then standard definition in terms of 
intcgro-differential operators. Geometrically, Guthrie's recursion operator for 
an equation £ is a Backlund autotransformation for the linearized equation V£ 
([22]). 

In geometrical terms, the linearization V£ can be introduced as the vertical 
vector bundle V£ — ► £ with respect to the projection £ — ► M on the base 
manifold. 

At the level of systems of PDE, the linearized system is 
^=0, l F i[U]=0, (18) 
where 

M£/] = £|^ (19) 

kj 1 

(cf. the Frechet derivative [28]), where U k are coordinates along the fibres of the 
projection V£ — ► £ and serve as additional dependent variables ('velocities' of 
the u k 's). We assume summation over all k, I such that the functions F l depend 
on Uj. 

Morphisms £ — > V£ that are sections of the bundle V£ — ► £ are in one- 
to-one correspondence with local symmetries of the equation E. Recall that 
nonlocal symmetries (more precisely, their shadows [15]) correspond to mor- 
phisms £ — ► V£ over £, where £ is a covering of the original equation. In full 
generality, Guthrie's definition includes such a covering. 

Let us denote by V£ — ► £ the pullback of the vertical bundle V£ — ► £ 
along the covering projection £ — > £. Then nonlocal symmetries correspond to 
morphisms £ — ► V£ that are sections of the projection V£ — ► £ . In coordi- 
nates, if the covering £ is determined by equations z 3 x = f J , z 3 y = g 3 , then its 
linearization V£ corresponds to the system 

F l =0, zi = f, z{ = g\ £ pl [U}=0. (20) 

Definition 22 ([11]) A recursion operator for the system (1) consisting of equa- 
tions F l = 0, I = 1, . . . , s, is given by the following data: 

(1) a g[„-valued zero-curvature representation a — A dx + B dy for £ ; 
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(2) two n- vector- valued functions A = (Ai), B Q = (Bl) on V£ linear on the 
fibres (i.e., linear in the variables Uj) and satisfying 



{D v - B)A a = (D x - A)B ; 



(21) 



(3) an s x n-matrix- valued function C on £; 

(4) an s- vector- valued function C D on V£ linear on the fibres (i.e., linear in 
the variables Iff). 

The following condition is supposed to hold: If U — (U k ) satisfies the linearized 
equation VE, then so does U' = L(U), where L(U) 1 = C)W^ + C D and W* , 
j = 1, . . . , n, are nonlocal variables of the covering 

W 3 X = A\ W* + Al U;; B;\V ■ //'. (22) 

The recursion operator defined by these data will be denoted as LK^ 1 . 

Once a is a ZCR and (21) holds, Equations (22) determine a covering; see [11, 
Eq. (3.2)]. 

Recursion operators exhibit the following form of gauge invariance: If S is a 
function on E with values in GL(n), then the data 



A' = A s = DxSS- 1 + SAS~\ A' = S*,4 , 
B' = B s = DySS- 1 + SBS-\ B' = SB , 
C = CS \ C = C 



(23) 



(we assume matrix operations) determine the same recursion operator as a map- 
ping U^U'. 

Remark 23 One can put the definitions in a more compact form. Let us con- 
sider (1 + n) x (1 + n) matrices 



/o 


^ 




(° 


\ 






, B = 






A 


A 




B 


B 


V 


J 




\ 


) 



(24) 



Then 

a = Adx + B dy 

is a ZCR for V£ ; this follows easily from formulas (21). Moreover, let us intro- 
duce the s x (1 + n)-matrix 



C= Co 



c 



Then the above formulas (23) of gauge invariance assume the compact form 
A' = A S = Djs- 1 + SAS~\ 

B' = B S = DySS' 1 + SBS~\ (25) 

C 1 = cs- 1 
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where 



A 


o \ 





s 


V 


J 



It is even possible to define a generalized recursion operator of the system (1) 
as consisting of a gljy-valued zero-curvature representation^ = Adx + B dy 
for V£ along with an s x N- matrix- valued function C on V£, subject to the 
following condition: If 



W{ = B\W\ 



(26) 



then U n = C l jW^ satisfies the linearized equation V£ . 
For A, B given by formulas (24), the correspondence between W and W is 



/ 7 \ 



W = 



V J 



where 7 satisfies D x j = D y -f — 0. With S being an arbitrary matrix, formulas 
(25) define a generalized gauge invariance of generalized recursion operators. 

Coverings (22) with a = are associated with conservation laws, since for 
them Eq. (21) reads D y A = D X B . Examples are provided by recursion oper- 
ators that can be written in the traditional integro-differential form ([27]) 



U" 



8 = 



Upon the obvious identification D I x U k = Uf and introduction of nonlocal vari- 
ables Wi — D~ 1 p 3 k U k , the Guthrie form of this operator is 

wi=tfui 

U n = C\Wi +R l k I U?, 

where pj!Uj dx + q^Uj dy is a conservation law of the linearized equation V£ 
(typically a linearized conservation law of the equation £; [22]). 

Example 24 The Lenard recursion operator D xx + Au + 2u x D x 1 for the KdV 
equation u t = u xxx + 6uu x has the following Guthrie form (with £ = £ and 
V£ = V£): 



W X = U, 

w t = u xx 
u' - u xx 



6uU, 
AuU- 



(27) 



2u x W. 
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Indeed, if U satisfies the linearized equation V£ , i.e., 



U t = U xxx + 6uU x + 6u x U, (28) 

then so does U' (for the same u). 

Here W is a potential of the conservation law U dx + (U xx + 6uU) dt of V£, 
which is a linearization of the conservation law udx + (u xx + 3u 2 ) dt of £. 



6.1 Inversion of recursion operators 

A recursion operator is said to be invertible if the morphism L of Definition 22 is 
a covering. The^recursion operator LK^ 1 is then simply a pair of linear coverings 
K, L : TZ — ► V£, its inverse being the recursion operator KL~ X . Noninvertible 
recursion operators do exist, see Remark 27(2). 

One immediately sees that a recursion operator and its inverse are built upon 
one and the same covering £ . In practice usually £ = £ (hence the covering £ 
is almost obsolete in Definition 22); however, one can simplify the ZCR a with 
the aid of a suitably chosen covering. Namely, given a recursion operator 

TZ 

V£ V£ 

associated with a ZCR a, the obvious pullback along a covering p : £ — ► £ 
yields a recursion operator 

p*TZ 

V ' K / \ P ' L 
V£=p*V£ P *V£ = V£, 

which is associated with the pullback p*a. 

For instance, let £ — > £ be the trivializing covering for a. Then, after suitable 
transformation (23), we have p*a = 0, whence the recursion operator becomes 
intcgro-diffcrcntial of first order in D _1 . Hence a possible way of conversion of 
recursion operators from Guthrie's to Olver's form, mentioned in the Introduc- 
tion. This approach was used in the work by Guthrie and Hickman [12] who, 
by using formal power series in the spectral parameter A, were able to describe 
large algebras of nonlocal symmetries of the KdV equation resulting from iter- 
ated application of the inverse recursion operator. 

Alternatively, £ — ► £ can be a covering such that p*a is strictly lower trian- 
gular (belongs to t^). Then the covering (22) will be abelian by a similar argu- 
ment as in Proposition 5 and the recursion operator will be integro-differcntial 
of order < s in £) _1 . 

Let us now turn back to recursion operators LK~ X with a general covering a. 
One usually observes that for systems £ integrable in the sense of soliton theory 
the covering if is of a very special form, which is described in the following 
proposition: 

Proposition 25 Let a = Adx + B dy be a Q-valued ZCR of equation £. Then 
the trivial vector bundle q x V£ — ► V£ carries a covering structure determined 
by the condition that an arbitrary element W of the Lie algebra g be subject to 
equations 

W X = [A,W}+£ A [U], W y = [B,W]+£ B [U}. (29) 
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Otherwise said, the associated ZCR a coincides with the adjoint representa- 
tion of the ZCR a, while A = £ A [U], B = l B [U\. 

Proof The validity of formulas (21) follows from the fact that A i — > ^[C^] is 
a differentiation. 

Taking account of the last proposition, we arrive at the following construc- 
tion, which converts a recursion operator from Guthrie's form to Olver's form 
provided the covering K is of the type (29). 

Construction 26 Step 1. Construct the generalized Riccati covering (Defi- 
nition 19) £' over £ such that a' := a H is lower triangular, where H is the 
matrix (11). 

Step 2. Let a' u , b' u be the diagonal entries of the lower triangular matrices 
A H , B H , respectively. Then a ' u dx + b' u dy are conservation laws; if they are 
nontrivial, then construct the abclian covering £" over £' with the corresponding 
potentials Zj. 

Step 3. Compute S = ZH, where Z is the diagonal matrix diag(e _Zi ). Obvi- 
ously, a" :— a s is then strictly lower triangular, and so is its adjoint represen- 
tation 

i? =a S , 

where S is the image of S in the adjoint representation of the group G. The re- 
part of the resulting recursion operator given by formulas (23) will be expressible 
in terms of inverse total derivatives D~* . 

Step 4 (optional). Let us consider the compact form (24) of the recursion 
operator, which now takes values in the algebra t„+i of strictly lower triangular 
matrices of dimension n + 1. Choosing appropriately a lower triangular gauge 
matrix S with units on the diagonal, one can, in principle, further simplify the 
formulas. 

If omitting Step 2, the recursion operator will be expressible in terms of in- 
verses (D x — a-j) -1 . 

Remark 27 (1) Let R be a conventional recursion operator of an integrable 
system, let id denote the identity map. As a rule, the inverse recursion operator 
(i? + Aid) _1 in the Guthrie form includes a A-dependent ZCR a. The parameter 
A can be usually identified with the spectral parameter of the standard ZCR of 
the system. 

(2) Let us recall that the formulas (29) can serve as a starting point of a 
method to find the inverse recursion operator of an integrable system without 
finding the conventional recursion operator first. One simply computes all mor- 
phisms 1Z — ► V£, where 1Z is the covering determined by (29). Recently the 
procedure has been applied to the stationary Nizhnik-Veselov-Novikov equa- 
tion, see [24]. Remarkably enough, the so obtained recursion operator turned 
out to be noninvertible for the zero value of the spectral parameter A. Two 
examples of such computation can be found below. 

7 Examples 

Example 28 Continuing Example 24, let us invert the Lenard operator. The 
result is, of course, well known (Guthrie and Hickman [12], Lou [20, 21]). 
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Instead of tedious inversion of the operator given by formulas (27) and (28), 
we compute it from scratch. We start with the standard s^-valued ZCR 

«=f ? nWf I* Uxx + 2U )dy 





1 0/ \— 2m -«j 

of the KdV equation with the spectral parameter set to zero. Using (22) and (29) 
with 5 [2 parametrized as 

Q P 
R -Q 

we get the following formulas for the covering K: 

Q 
'P 

Q J = \ 2u u xx + 2u 2 

K Rj t \ -4m -2u x 

Here U denotes a symmetry of the KdV equation, i.e., satisfies the linearized 
KdV equation (28). Then one easily finds that U' = Q satisfies the same lin- 
earized KdV equation (28) as well, i.e., yields a recursion operator for the KdV 
equation. It is a matter of routine to check that this operator is the inverse of 
the Lenard operator. Moreover, it follows that K : 1Z — ► V£, originally given 
by U — U' xx + AuU 1 + 2u x W, constitutes a three-dimensional covering (with 
nonlocal variables U, U x and W) . 

According to Construction 26, to express the inverse recursion operator in 
terms of D~ x , all we need is to make the ZCR a strictly lower triangular. As 
the first step we build up a covering £' — ► £ with the quadratic pseudopotential 
h = fin defined by Eq. (14), i.e., 

h x = —h 2 — m, 

h t = —2uh 2 + 2u x h — u xx — 2m 2 . 
Then, using the gauge matrix 

we get the lower triangular ZCR 

\-l hj \ -2m -u x + 2uh) y 

with —h, h on the diagonal. As the second step, we construct the abelian covering 
£" — ► £' with the potential z satisfying 

z x = —h, z y = u x — 2uh. 

The gauge matrix 
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then gives the strictly lower triangular ZCR 

a" = a ZH =( ° 2z Jote+l o ° 22 ° n )dy. 



2 Z ; "V-2c 2 ^ 

In the third step, we combine the above gauge matrices into one and compute 
its adjoint representation: 

, _ 7 , _ 7 . /c- 2z -2he- 2z -h 2 e- 2z ^ 

V 7 \ <?' 

Acting by S on our operator, we get 
'P\ 

Q 

K RJ 
'P\ 

Q 

2e~ 2z hU' x + (2h 2 + 4u)c- 2z U'^ 
U' x - 2hU' 
-2c 2z U' 

U = Q- hc- 2z R. 

Rewritting the x-part in terms of inverse total derivatives D~ x , we get P - 
D-\er 2 *U), Q = D-^e^'P), R = -2L>- 1 ( e - 2z Q), hence 

U = D-^D^e-^U - he- 2z D- 1 e- 2z D- 1 e 2z D- 1 e- 2z U. 

This is the well-known result [12, 20, 21], since U' = Q-hc- 2z R = -\D x (R/h 2 
and z xx = z 2 +u. 
The optional fourth step does not bring any significant improvement. 

Example 29 Let us consider the Tzitzeica equation [38] 




-2m 



later rediscovered as a member of the Zhiber-Shabat classification [45] . Its ZCR 
l-u x \\ 

a = \ A u x | dx+ ( e u /X]dy (30) 

A Oy 

as well as the Backlund transformation were essentially found by Tzitzeica him- 
self. 

One could invert the known recursion operator [36], but it is easier to com- 
pute the inverse recursion operator directly by the procedure outlined in Re- 
mark 27(2). Namely, we consider the eight-dimensional covering (29), where 
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A,B,A Q and B Q are found from the formula (30) to be 
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-A 
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-2A 





-u x 
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u x 
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u x 
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-u x J 






( 








— e 


U /A e" 


- 2u /x 













° \ 
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-e"/A 















e" 2 "/A 


























-e"/A 


e"/A 
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-2uj x 













c u /A 




-e"/A 













-2e u /A 













2e u /A 













e u /A 
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e"/A 


















-c- 2u /X 


o J 


(- 


-u x \ 
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\ 











A a = 







u x 




V o / 



B = 



V 



2e~ 2u U/X 




e u U/X 
e u U/X 




W being a column (Wn, W12, W13, W21, W22, W23, W31, Ws2) T of pseudopoten- 
tials. One easily finds that Wn — W22 is a symmetry of the Tzitzeica equation 
if so is U. We have obtained the 'inverse' recursion operator of the Tzitzeica 
equation in the Guthrie form. 

Let us express it in terms of D^ 1 . As the first step we introduce pseudopo- 
tentials p, q, r satisfying 



p x = Xp 2 - 2pu x - Ag, 

q x = Xpq - qu x - A, 

r x = -Xpr + Xq + Xr 2 + u x r, 



Py 



1 



-pq 



e A 



q y = y(? 2 -p), 

r y = —{-pr 2 + qr - 1). 



to make the ZCR (30) lower triangular by providing a solution to Equations (17). 
Indeed, acting on a by the gauge matrix 




22 



we get 



H 



-u x + Xp 

( i " = | A u x — \p + \r I fix 
A -Ar 

>e u q/\ 
e"r/A -e u pr/X | dy. 

,e"/A -e>/A e u (pr - q)/\ 

In the second step we remove the diagonal. To this end we introduce pseudopo- 
tcntials s, t by 



s x = ~u x + Xp, 
t x — u x - Xp + Ar, 



-pr. 



Acting on a H by the gauge matrix 

we finally get 

„zh 



-» 






















e s+t 




Xe s+2t 



+ 





a u+s— t 



„u+2s+t 



r/X 
/A 





-e""°'"V/A Oy 

Denoting S = ZH, we compute the adjoint representation S to be 



2s-t 


-e- 2s -*r 


c- 2s - l p e" 2s 


-\pr-2q) 


-e- 2s ~\pr + q) 





e -s+t 





-e~ s+t p 


e- s+t p 








e -s-2t 


c s-2t r 


-2c- s - 2t r 
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' 2s ^p{pr - q) 


~ e - 2s -y 


e'^-tqipr - q)\ 






-e- s+t p 2 


e~ s+t q 


-e' s+t pq 




e 


- s - 2t (pr-q) 


-r 2 e- s - 2t 


e~ s ~ 2t (pr - q)r 






P 





q 






-P 
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—pr 
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e s -*r 
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-c s+2t p 












c 2s+t ) 
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Acting by S on the above recursion operator we get 






























\c s+2t 
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-Ae s "* 
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-Xc s+2t 




























Ae 5 "' 


-Ac s - ( 






















Ae*+ 2 


2Xc s+2t 











V 

















Ac s+2t 


-Ae s -* 


0/ 



and 

/e^-'Mpr + q ) Ux \ 
2e- s+t pU x 
-e- s ~ 2t rU x 






V / 

(we omit the matrices B s and B Q ). 

Thus, the inverse recursion operator for the Tzitzeica equation in terms of 
D- 1 is 

V = W 2 i - W 22 - 2c- s+t pW2 3 + c- s ~ 2t rW 31 + e- 2s - t (2pr - q)W 32 , 



where 






W u 


= D~ 


1 [e- 2s - t (-2pr + (7 )C/J, 


w 12 


= D~ 


1 [2e~ s+t pU x ~e s+2t XW 11 ], 


w 13 


= D~ 




W22 


= D~ 


1 [U X + e s -*AWi2 - e s+2t A^i3 


W21 


= D~ 


1 [-u x - c s - t \w 12 ], 


w 31 


= D~ 


1 [e s+2t A(VK 2 i + 2VK 22 )], 


W23 


= D~ 


1 [-e s - t \{~W 2 i + W22)}, 


W 32 


= D~ 


W+^-e^VML 
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